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ABSTRACT 

This paper presents the results from the Jodrell Bank - lAC two-element 33 GHz 
interferometer operated with an element separation of 32.9 wavelengths and hence 
sensitive to 1° -scale structure on the sky. The level of CMB fluctuations, assuming a 
flat CMB spatial power spectrum over the range of multipoles I — 208 ± 18, was found 
using a likelihood analysis to be AT^ — 63l'lg /iK at the 68% confidence limit, after the 
subtraction of the contribution of monitored point sources. Other possible foreground 
contributions have been assessed and are expected to have negligible impact on this 
result. 

Key words: cosmology: cosmic microwave background - cosmology: observations - 
large-scale structure of the Universe - instrumentation: interferometers. 



1 INTRODUCTION 

Observations of the angular power spectrum of the CMB 
temperature fluctuations are a powerful probe of the fun- 
damental parameters of our universe. The amplitude and 
spatial distribution of these fluctuations can discriminate 
between competing cosmological models. Most inflationary 
models predict more power on scales of 0? 2 - 2°, in the form 
of a series of peaks. These are due to acoustic oscillations 
in the photon-baryon fluid, which are frozen into the CMB 
at recombination, with the peaks corresponding to regions 
of maximum compression and troughs regions of maximum 
rarefaction. Hence, the position of the flrst acoustic peak is 
a strong test for the geometry of the universe, since it corre- 
sponds to a fixed physical scale at the time of recombination 
projected onto the sky. 

The previous result from the Jodrell Bank - lAC 
33 GHz interferometer of ATf = 43^12 reported in 
Dicker et al. (1999), corresponds to an angular spherical 
harmonic I ^ 110, equivalent to ~ 2° structure. To investi- 
gate smaller angular scales the baseline was doubled; in this 
paper we analyze the data from this wide spacing configu- 
ration which corresponds to an angular spherical harmonic 
I ~ 210. The data presented here were taken at the Teide 
Observatory, Tenerife, between 27 May 1998 and 9 March 
1999. The paper is organised as follows. The instrumental 
configuration is summarised in Section |^ a full description 
can be found in Melhuish et al. (1999). The basic data pro- 



cessing is outlined in Section for a more complete discus- 
sion see Dicker et al. (1999). The calibration method is also 
discussed in Section g and the data analysis in Section ^ 
A derivation of the fluctuation amplitude, after an estimate 
of the contribution of possible foregrounds, is given in Sec- 
tion ^ 



2 THE 33 GHZ INTERFEROMETER 

The interferometer consists of two horn-reflector antennas 
positioned to form a single E-W baseline, which has two 
possible lengths depending on the separation of the horns. 
The narrow spacing configuration has a baseline of 152 mm 
while in its wide spacing configuration the horns are 304 mm 
apart. For the observations presented here, the baseline was 
304 mm. Observations were made at a fixed declination of 
Dec=-|-41° , using the rotation of the Earth to "scan" 24*" in 
RA each day. This "scan" runs through some of the lowest 
background levels of synchrotron, dust and free- free emis- 
sion. The horn polarization is horizontal - parallel with the 
scan direction. There are two data outputs representing the 
cosine and the sine components of the complex interferome- 
ter visibility. The operating bandwidth covers 31 - 34 GHz, 
near a local minimum in the atmospheric emission spectrum; 
the antenna spacing corresponds to 32.9 wavelengths. The 
low level of precipitable water vapour, which is typically 
around 3 mm at Teide Observatory, permits the collection 
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Figure 1. The total data stack, in 0.5° -bins, at Dec +41° collected over the period 27 May 1998 to 9 March 1999 and showing the 
cosine and sine visibility data. The number of days' observations included at each RA is between 180 - 210 days due to the removal of 
the Sun during the year. The Galactic plane crossing at the anti-centre is clearly visible at 21*^ together with Cyg A at 20*^30™ and 3C84 
at 3^20™ 



of high quahty data. Only 16 per cent of the data have been 
rejected, due to bad weather and the daily Sun transit. 

The measured response of the interferometer is well ap- 
proximated by a Gaussian with sigmas of ora= 2? 25 ± 0?03 
(in RA) and ODec= 1°00 ± 0?02 (in Dec), modulated by 
fringes with a period of / = 1?74 ± 0?02 in RA. This de- 
fines the range of sensitivity to the different multipoles i of 
the CMB power spectrum (Ce) in the range corresponding to 
a maximum sensitivity at £ = 208 (0?8) and half sensitivity 
at Ae = ±18. 

A known calibration signal (GAL) is periodically in- 
jected into the waveguide after the horns allowing a con- 
tinuous calibration and concomitant corrections for drifts in 
the system gain and phase offset. 



3 BASIC DATA PROCESSING AND 
CALIBRATION 

The first step in the analysis is the removal of any variable 
baseline offsets from the data and the correction of a small 
departure from quadrature between the cosine and sine data. 
The data are calibrated relative to the GAL signal and re- 
binned into 2 - minute bins to ensure alignment in RA be- 
tween successive scans. The data affected by the Sun and 
bad weather are removed and individual scans are weighted, 
with respect to their RMS error, to form a "stack" . The to- 
tal stack of all the data used for this analysis is shown in 



Fig. 0. The number of observations at each RA varies from 
210 to 180 days due the removal of ±0.5'^ about the Sun 
transit each day. 

The data are calibrated relative to GAL, although GAL 
itself needs to be calibrated by an astronomical source. The 
large size of the primary beam results in a reduced sensitivity 
to point sources and many days of observation are required 
to achieve a signal - to - noise ratio sufficient for calibration 
purposes. Gonsequently, the Moon is used as the primary 
calibrator as the power received from a single Moon transit 
is large enough to give signal - to - noise ratios of '-^ 6000. 

The Moon was modelled as a uniform disk of radius 
T'Moon and a 33 GHz brightness temperature, Tb given by: 



Tb = 202 -I- 27cos( 



K 



(1) 



where (f) is the phase of the Moon (measured from full Moon) 
and e = 41° is a phase of fset caused by the finite th ermal 
conductivity of the Moon (Gorenstein & Smoot 1981). It is 



sufficient to model the Moon as a uniform disk when correct- 
ing for its partial resolution by the interferometer. The ef- 
fect of temperature variations across the Moon are negligible 
when compared to the error, 5.5 per cent, of the Gorenstein 
and Smoot model for the Moon's brightness temperature. 
The expected antenna temperature, Te, can then be found 
by integrating over the disk of the Moon, multiplied by the 
normalized interferometer beam function: 
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where x = rcostp and y — rsin?/;. ctra and obcc are the 
RA and Dec beam sigmas (dispersion) and / is the fringe 
spacing. 

Regular observations of the Moon were made; for each 
observation, equation ^ was evaluated numericaUy and an 
amphtude for CAL found such that the amphtude of the 
Moon in the processed data was equal to the predicted value. 
Using 27 observations of the Moon, an average amplitude for 
CAL of 14.7 ± 0.8 K was found. The error consists of a 1.4 
per cent error in the meas urements and an estimated 5.5 per 
cent in the moon model (Goronstein & Smoot 1981). Fig. ^ 



shows our observations and how the measured brightness 
temperatures of the Moon change with phase. 

The moon model used in this paper differs from the 
model given in Dicker et al. (1999). This data together with 
the addition data taken in the narrow spacing, section ^ 
will be the subject of a forthcoming paper, using the Moon 
model given in this paper. 



4 LIKELIHOOD ANALYSIS 
4.1 Theory 

The temperature anisotropies of the CMB fluctuations are 
described by a two dimensional random field on the sky, the 
properties of which can be determined from the two point 
correlation function C'^^^ {9ij), where 6ij is the angular sep- 
aration of the two points. 
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Figure 2. Measured value of the Moon's brightness temperature 
as a function of phase. Phase is measured from the full Moon. The 
solid line is the predication of the model given by Gorenstein & 
Smoot (1981): Tb = 202+27 cos((/)-e) in Kelvin. Each observation 
has been calibrated using CAL = 14.7 K. 
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which can be expanded in terms of spherical harmonics as; 
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In our analysis the form of the likelihood function is given 
by: 

1 



L oc 
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-exp 
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(5) 



where D is the data set and C is the covariance matrix, 
which represents the model of the CMB sky modulated by 
our observing strategy. The covariance matrix is composed 
of two terms, C = S + N where S is the signal and TV is 
the noise correlation matrix. The signal is the convolution of 
the two point correlation function and the auto - correlation 
function of the primary beam function of the interferometer. 



Si' 



obeam / n 



(6) 



U sing the band 
^l{i+l)Ci/2-K I 
range of I covered by the window function 



power approximation where ATf = 
is assumed to be constant across the 



Figure 3. The window function of the interferometer in its wide - 
spacing configuration (304 mm ^ 32.9 wavelengths)which is well 
fit by Wm = 0.037 X exp ( - i£^H2|i£^ ) 
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where Imax is the limit of the summation. 

The sensitivity at an individual value of I is given by the 
window function of the interferometer in Fig. ^. The window 
function was computed using the method of Muciaccia, Na- 
toli & Vittorio (1997) to decompose the beam into spherical 
harmonics. The resulting function can be modelled by: 



WfXQ) = 0.0377 X exp 



(^- 208)(^- 209) 
694 



(8) 



Wi{Q) falls to half power at ^ = 208 ± 18. 



* There is a typographical error in Dicker et al. (1999), in the 
denominator the Stt should be 2it as written here. 



The covariance matrix can now be calculated, by nu- 
merically evaluating equations (pi) and (|7|). Equation (fTl) is 
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Table 1. Sources within a 4°-wide Dec strip centred on +41°. 



Name RA2000 Dec2000 Mean Flux (Jy) 



3C 84 


031^ 


19™ 


48= 


+41° 30' 


42" 


11.8 


DA 193 


05*' 


55™ 


31= 


+39° 48' 


49" 




4C 39.25 


09^ 


27™ 


03= 


+39° 02' 


21" 


9.5 


3C 345 


161^ 


42™ 


59= 


+39° 48' 


37" 


9.0 


BL Lac 


22^ 


02™ 


43= 


+42° 16' 


40" 


3.7 



evaluated to Imax ~ 300, where the sensitivity of the inter- 
ferometer, W^^(O), is neghgible. The resuh of equation (|^, 
C'^^^{9ij), is used to calculated the required value of Sij 
using equation (^) . The covariance of two points with an 
angular separation of Oij is then given by dj = Sij + Nij. 

4.2 The results and the Galactic Cut 

Any analysis should take account of likely Galactic emission. 
At Dec +41° the ranges 2l'^48™- 3*' 48™ RA and 6''12™" 
19''30°' RA are at Galactic latitude b > 10° . To find re- 
gions free of significant Galactic emission, 48 intervals of 
5 hours in RA were analysed as described above, stepping 
every 0.5 hours. It was found that the above cut was in gen- 
eral agreement with the higher latitude results, except the 
region 2''48™- 3''48™ which was excluded; this is due the 
Galactic plane emission and diffi cult ies in subtracting the 
contribution of 3C84, see Section 5.1. 

The likelihood analysis of the ranges 2l''48™- 2''48™ 
RA and 6''12™- 19''30™ RA gives AT = 78.5ti2:l ^iK 
and AT = 69.5li[2 q ^oi" cosine channel and the 
sine channel respectively. The likelihood analysis using both 
channels simultaneously gives AT = 70.0^g 5 at the 68 
percent confidence level. 



5 THE EFFECT OF FOREGROUNDS ON THE 
RESULTS 

5.1 Point Sources 

The 5 strongest sources with 5(33 GHz) > 2 Jy within a 4° 
strip centred on Dec +41° , listed in table 0, are routinely 
monitored by the University of Michigan at 4.8, 8.0 and 
14.5 GHz and in the Metsahovi programme at 22.0 and 37.0 
GHz. Using these data over the period of our observations, 
it was possible to assess their mean fiux densities at 33 GHz. 
These were then convolved with the two-dimensional inter- 
ferometer beam pattern centred on Dec +41° and converted 
to antenna temperatures using the factor 6.90 /^KJy~^; in 
this form these sources may be subtracted from the data. 

The data ranges 21'^ 48™- 2'M8™ RA and 6''12™- 
19''30™ RA were analyzed together, subtracting the point 
sources as discussed above. Each channel was analyzed in- 
dependently and then combined for a joint analysis. For 
the cosine channel AT = 69.5tii j ^K, for the sine chan- 
nel AT — 62. 515^1 5 and combining both channels gives 
AT = 64.0j;6:o mK at the 68 per cent confidence level and 
AT = QA.Otlt l mK at the 95 per cent confidence level. The 
likehood curve of this anaylsis is shown in Figure ^. 




100 



AT(/xK) 

Figure 4. The likelihood curve of the joint analysis of both chan- 
nels of the point source-subtracted data. AT = 64.0^g g fiK at the 
68 per cent confidence level and AT = 64.0^j^3 Q fiK at the 95 per 
cent confidence level. 



The contribution of unresolved point sources was esti- 
mated according to the results of Franceschini et al. (1989), 
at 33 GHz resolution of 0.8° this is expected to be ~ ll^K, 
which adds in quadrature to the CMB signal. The con- 
tribution of unresolved sources then accounts for approxi- 
mately 1 fjK of the total signal and accordingly the best es- 
timate of the intrinsic CMB fluctuation amplitude is AT = 
63.0lg;U ^iK at ^ = 208. 

5.2 Spinning Dust 

de Oliveira-Costa et al. (1998) estimated the contribution of 
spinning dust in a 19 GHz map of 3° resolution by correlat- 
ing it with the DIRBE sky maps, finding AT ~ 66 ± 22^K. 
Using the IRAS 100 fiK map, Gautier et al. (1992) investi- 
gated the spatial index of the dust, finding on scales between 
8° and 4' that AT cx £~^^^. This, combined with the ex- 
pected spectral index of the spinning dust of —3.3 < /3spin < 
—4 (de Oliveira-Costa et al. 1998); allowed the estimation of 
the contribution of spinning dust at 33 GHz and 0.8° reso- 
lution. The expected signal in our data due to spinning dust 
was found to be AT'^"=* « 1.5 ± 0.5/iK. This again adds 
in quadrature to the total signal, therefore the contribution 
from spinning dust is expected to be negligible. 



5.3 Diffuse Galactic Emission 

An estimate of the amplitude of the diffuse Galactic compo- 
nent in our data can be computed using the results obtained 
in the same region of the sky by the Tenerife CMB experi- 
ments (Gutierrez et al. 2000). At 10.4 GHz and on angular 
scales centred on ^ = 20 the maximum Galactic component 
was estimated to be < 28 ^K. Assuming that this contribu- 
tion is entirely due to free-free emission (/3 — —2.1) and a 
conservative Galactic spatial power spectrum of £~^'^, the 
predicted maximum Galactic contamination in the data pre- 
sented here is 0.8 fiK, less than 2 per cent of our measured 
value. Any such contribution would add in quadrature to 
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Figure 5. Values of AT published in the past year, as a function of I. Our present result of AT^ = 63^g /xK is shown by the heavy 
line at £ ~ 200. The result of Dicker et al. (1999) is shown by the heavy line at £ ~ 100. The dotted line represents the model given by 
Ob = 0.05, f2cDM = 0.40, Q,x = 0.55 and Hq = 70kms~^ Mpc~^; this is shown for illustrative purposes and does not represent a fit to 
the data. 



Table 2. Values of AT around £ ~ 200 published in the past year. 



Experiment Freq. (GHz) AT (pK) Multipole, £ Reference 



Int33 


33 


631^ 


208+^1 


This Paper 


MSAMl 


156 




201+^0 


Wilson et al. (2000) 


Viper 


40 


66_i7.2 


23711^1 


Peterson et al. (2000) 


Boomerang 


150 


79 + 10 
'^-10 


204^2^ 


Mauskopf et al. (1999) 


Python V 


41 




199+^5 


Coble et al. (1999) 


TOC098 


144 


83+1 


226+37 


Miller et al. (1999) 


TOC097 


35 




199+_ll 


Torbet et al. (1999) 



that from the CMB, and so is insignificant. The true make- 
up of the 10.4 - GHz Galactic foreground emission will have 
a steeper average spectral index since synchrotron radiation 
with /3 ~ — 3 will contribute to the measured value, there- 
fore the contribution to our result will be even lower than 
stated. 



6 CONCLUSIONS 

In this paper we describe the results taken with the Jodrell 
Bank - lAC 33 GHz interferometer in its wide spacing con- 



figuration, which is sensitive to structure aX I — 208 ± 18. 



In the final result of AT = 63.01g Q /iK possible foreground 
contributors have been considered, the most significant of 
which are point sources. In Section 
strong, monitored sources are removed from the data. 



the contribution of 
To 

allow for the contribution of unresolved sources the lower 
limit on AT has been increased. 



and 



5.3, the 



In Sections 5.2 

contributions from dust and diffuse Galactic emission are 
found to be negligible. We believe our result represents the 
intrinsic CMB fluctuation signal at £ — 208±18. The quoted 
error in our result is dominated by sample variance, result- 
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ing from the finite number of beam areas observed. 

Table ^ and Fig. ^ sliow our result alongside others pub- 
lished in the last year from experiments covering similar an- 
gular scales. These experiments have been made at a range 
of different frequencies and regions of the sky. Our result is 
in good agreement with the published data around £ ~ 200. 
The results at ^ ~ 200 appear to be converging on a value 
of AT = 60 — 70 /iK. However, over the wider-^ range of 
50 - 300 discrepancies appear in the data sets significantly 
greater than the quoted errors. There appears to be evidence 
for unknown systematic effects and possible foreground con- 
tamination remaining in the data sets. 

Our interferometer results show a rise in the amplitude 
of the power spectrum between £ ~ 100 and £ ~ 200. This 
is intrinsic to the CMB since the parameters of the inter- 
ferometer system, particularly the calibration, remain the 
same except for the spacing. The data in Fig. BI despite the 
discrepancies referred to above, are strongly indicative of a 
peak in the power spectrum at £ ~ 200. 

The interferometer is currently in its narrow spacing 
configuration {£ ^ 110), observing declinations spaced by 
1?2 from Dec -1-37.4 to -f43.4; these data will significantly 
reduce the sample variance of the result published by Dicker 
et al. (1999) to the order of 5 per cent. 



Recontre de Moriond, Microwave Background Aniso tropies. 
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Note added in proof 

Our rcsuh of AT = GS.OtH )j,K at £ = 208 is competitive 
with the latest published results from the Boomerang ex- 
periment, which found an amplit ude for the first peak o f 
AT = 69 ± 8 at £ = 197 ± 6, ( |de Bernardis et al. 2o6c| ) 
and the Maxima experiment, w hich found a peak at £ ~ 220 
of amplitude AT = 78 ± 6 /xK, ([Hanany et al. 2000|). 
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